Effective free energy for pinned membranes 
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We consider membranes adhered through specific receptor-ligand bonds. Thermal undulations of 
the membrane induce effective interactions between adhesion sites. We derive an upper bound to 
the free energy that is independent of interaction details. To lowest order in a systematic expansion 
we obtain two-body interactions which allow to map the free energy onto a lattice gas with constant 
density. The induced interactions alone are not strong enough to lead to a condensation of individual 
adhesion sites. A measure of the thermal roughness is shown to depend on the inverse square root 
of the density of adhesion sites, which is in good agreement with previous computer simulations. 
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Living cells interact with their environment through 
the cell membrane, a bilayer of phospholipids consisting 
of hydrophobic tails and hydrophilic heads. Embedded 
integral proteins, which make up a large portion of the 
membrane, are involved in a range of processes such as 
ion passage and cell signaling. Through the formation 
of receptor-ligand bonds certain proteins such as cad- 
herins, selectins, and integrins are involved in cell adhe- 
sion, where one cell binds to another cell, to an extracel- 
lular matrix, or to a surface [T]. Extensive studies have 
been conducted on cell mimetic model systems using a 
variety of techniques [SHS] . 

The proteins embedded into the membrane can move 
rather freely together with the lipids. The lateral dynam- 
ics and diffusion properties of embedded integral proteins 
have been studied both theoretically [HI E] and experi- 
mentally [5j. Protein dynamics is influenced by the cy- 
toskeleton, geometric constraints, and the thermal un- 
dulations of the membrane itself. Adhesion sites (i.e., 
single receptor-ligand bonds) in cells are not randomly 
distributed but form focal adhesions which are several 
square microns in area 0. Since the bending of the 
membrane requires energy one would imagine that such a 
clustering of adhesion sites is preferential. However, com- 
puter simulations of confined membranes [TUHO] indicate 
that fluctuation-induced interactions between adhesion 
sites are short-ranged and insufficiently strong to lead to 
thermodynamically stable clusters, quite in contrast to 
earlier predictions [T31 \T5\ . The aim of this Rapid Com- 
munication is to derive an effective and accurate expres- 
sion for the interaction free energy of pinned membranes 
through a variational approach. Membrane-induced in- 
teractions between adhesion sites are indeed found to be 
rather weak, confirming that active processes and inter- 
protein interactions are primarily responsible for the for- 
mation of adhesion domains. 

We consider a tension-less, almost flat membrane mov- 
ing close to a substrate in a fluid at temperature T. The 
membrane height profile h{r) is given in the Monge rep- 
resentation with r — {x, y). The projected area is and 
for simplicity we employ periodic boundary conditions. 
In this geometry the membrane can be thought of as a 
two dimensional sheet, where the thermal fluctuations 



are governed by the Helfrich energy [TB]. We focus on 
the bending energy 

no[hir)]= J d'r'^[V'h{r)r (1) 

as the dominant contribution with bending rigidity k. 
In the presence of N adhesion centers at locations {r^} 
pinning the membrane to the substrate the partition sum 
reads 

ZnHv,}) ^ / [dh{r)] e-««o+«.) J] eS{h{v,)) (2) 
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with f3 = (fce^)^^. The path integral sums over all ac- 
cessible membrane conflgurations h{v) ^ 0. The (small) 
length £ is related to the flexibility of bonds and quanti- 
fles the fluctuations of the membrane at the adhesion 
sites (one should think of the (5-functions as limits of 
Gaussians). Van der Waals forces, the repulsion due to 
the formation of an electrostatic double-layer, and steric 




FIG. 1: Sketch of a membrane with instantaneous height pro- 
file h{r) moving above a substrate. The membrane is pinned 
down through adhesion centers with effective size a. a) In the 
presence of only a few centers the thermal roughness away 
from the adhesion centers is ^x,d. b) For a large number of 
adhesion centers the membrane is pulled closer to the sub- 
strate. This reduces the thermal roughness < C±,d due to 
interactions of the membrane with the substrate. 
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forces give rise to an interaction energy 7ii[h{r)] between 
substrate and membrane, see Fig. [T] While explicit ex- 
pressions can be found for these forces [T7], here we pur- 
sue a more general strategy that does not depend on 
specific details. 

We expand the height profile h{r) into Fourier modes, 



where h* 



q 
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„q — .»-q. In particular, ho is the height averaged 
over the area of the membrane. The Hclfrich Hamilto- 



nian ([!]) becomes 



Ho — 



qeQ 



where we split = h'^ + \h'^ into real and imaginary 
part and q= |q|. We sum only over independent wave 
vectors, i.e., for a given q the set Q does not contain — q. 
Also, g = is excluded from the sum. 

Due to the pinning the membrane is pulled towards 
the substrate. The interactions of the membrane with 
the substrate lead to a 'confinement' of fiuctuations. For 
membranes placed between walls the interaction energy 
is often expanded into a quadratic function of the height. 
However, in principle we can do better through introduc- 
ing a quadratic reference energy Hj. that acts on every 
height mode independently. Employing the Jensen in- 
equality we know that the free energy is bounded from 
above as ITSi 



J'n^Fm + (AH) = F, 



N 



(3) 



with Tn = — ^InZ^r. Here, AH = Hi — Hi and 
(• • •) denotes the average with respect to the reference 
system with partition sum Zjy and free energy Fj^ = 
— /3~^lnZjv. By minimizing the function Fn with re- 
spect to a set of free parameters we obtain an approxi- 
mation of the true free energy J^jv({ri}). In this work we 
choose an intermediate description 
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with free parameters {70, /i, 7}. The total energy then 
reads 

mo+Hi) = ^%{ho -hf + ^J2 %K + (5) 

qeQ 

with 70 = Pjo and 7, = 2L^/?(7 + nq'^). 

Before we calculate the reference partition sum Zjy let 
us assume we had performed the minimization of Fn 
in the presence of the adhesion centers with constraints 
h{ri) = 0. We hypothetically remove the constraints but 
keep the values for 70, h, and 7 we had obtained as a 



result of the minimization. We define the thermal rough- 
ness ^± through 

ei ^ {[h{r) - ho]\ - Ed'^'il')" = 4 E (6) 

7? 
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which is independent of 70 and h. The brackets (• • •)o 
denote the average involving the energy Eq. ([s]) in the 
absence of constraints. For the sake of simplicity and for 
sufficiently large L we evaluate the sum over wave vectors 
q through a two-dimensional integral over the half-plane 
with the well-known result 
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dqq^~^l\/VS^^U,d (7) 
27r/L Iq 
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in the limit i ^ ^y. Here, = (k/j)^^'^ is the lateral 
decay length over which membrane fluctuations are cor- 
related. The thermal roughness is bounded by its value 
^_L.d ~ L/ y^lGn^/SK for a free, detached membrane with 
7 = 0. Note that the thermal roughness obtained by 
averaging over the projected area including the adhesion 
sites is greater than 

We integrate over the membrane height fiuctuations 
in the presence of N adhesion centers at positions 
{^i} [HI I19| . The reference partition sum including the 
constraints reads 



Zn = {£/2nf / [dh{r)] e-^(«o+«.) / dAi • • - dA 
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with coeSicients 
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Cq = 2^AjCosq-rj, ee 2 ^ A; sinq • 
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We have replaced the (5-constraints [Eq. ([2|] by integrals 
over a set of auxiliary variables {A^}. We assume that 
the excluded volume constraint is represented by the in- 
teraction energy Hi, and therefore by H^. Performing 
the integrations we obtain 



^ Zo[detm/^2Arj-i/2^ Zo=l[^. 

qeQ^-? 



(9) 



where we have dropped the sub-extensive contribution 
due to the integration over ho- The symmetric N x N 
matrix m has components = m{\ri — rj|), where 



qeQ 
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(10) 



describes the membrane-mediated interactions. The 
spatial correlations between height fluctuations and 
therefore the effective lateral interactions between ad- 
hesion sites is short ranged with potential u{x) = 



3 



— (4/7r) keio(a;), where keio(a;) is a Kelvin function pU] , 
The effective interactions decay on the single length 
scale ^1 that also determines the thermal roughness 
[Eq. §]. 

For the calculation of [Eq. ([s])] we first determine 
the contribution 



lnZo= }_^ln7, = iV(8Cif), 



(11) 



which for small captures the loss of entropy due to a 
large confinement of fluctuations. To establish F/v as a 
rigorous upper bound to the free energy J-^ we impose 



/3(AH) s=;-/3(H,)o = -iV(4e||)' 



This bound restricts the accessible values of 70 and h. 
We assume that the energy (Hr)o required to constrain 
the fluctuations is the dominant contribution to (AH). 
The equal sign holds if the interaction energy ("Hi) is ex- 
actly balanced by {Hi) — ('Hr)o describing the energy re- 
quired to bend (and, in principle, stretch) the membrane 
around the adhesion centers. Putting together Eqs. 
(12) we see that the resulting upper bound -F/v(C||; {ri}) 
is a function of alone and independent of interaction 
details between membrane and substrate, which is our 
first main result. In the following we use the minimum 
F;v(Cy ) as an approximation to the free energy. We split 
Fjv = Nf + into an ideal part / and an excess free 
Introducing an effective size a of adhesion 



energy i-jf 
centers with area fraction 
reads 



Na^/L^ the ideal part 



/3/(eii) 



(«/Cii)^ 
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+ ln(e||/a)+ln(a/eo) (13) 



with ^0 = VSfJni. The excess free energy becomes 



(14) 



due to the membrane-induced interactions between adhe- 
sion sites. These interactions are encoded in the matrix 
X with elements 



X^i(C||;{r,}) 




Using i < S,_L ^ C_L.d we find the meaningful range 



^0 < Cii sS 



1 



.L ~ 0.13L 



(15) 



(16) 



for the decay length of the effective interactions. For 
large ^y we expect a crossover to free fiuctuations while 
below ^0 the membrane is effectively so stiff that a trans- 
lation of adhesion centers does not require to bend the 
membrane. Using Eq. (12) the temperature and bend- 



ing rigidity /3k only set this crossover length and do not 
influence the value Q that minimizes F^. 
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FIG. 2: (color online) Left panel: the mean thermal roughness 
(^x) vs. area fraction for the lattice gas (circles) and the 
non-interacting adhesion centers [Eq. (17 1, solid line]. For 

m 



(12) comparision the global thermal roughness reported in Ref. 



is show (squares) with fit ^± ~ 0.140"^^'^ -I- 0.035 (dashed 
line). Inset: probability distribution (solid line) of for the 
lattice gas {(f) — 0.3), which is well described by a Gaussian 
(dashed line). Right panel: the mean effective interaction 
strength {v). The standard deviation of v is smaller than 
symbols. 



We first consider non-interacting adhesion sites, i.e., 
the off-diagonal elements of x vanish and F^ — 0. Min- 
imizing /(Cii) leads to the result 



a = 



(17) 



Such a scaling of the thermal roughness with the inverse 
square root of the area fraction has been found in com- 
puter simulations [10]. Moreover, the prefactor ~ 0.14 
for ^± determined numerically in Ref. |10| compares fa- 
vorably with our 1/8, see also Fig. |2] Of course, the 
divergence for </>—>■ is only apparent since we expect a 
crossover to free fluctuations with ~ ^j^.d if the mem- 
brane is pinned by only a few adhesion centers. 

In general ^ = ^{{'ri}) will depend on the geometry 
of adhesion sites and, therefore, is a fluctuating quantity 
for mobile sites, i.e., ligands are able to move in the sub- 
strate (e.g., through using a supported lipid bilayer as 
substrate). The effective interactions between adhesion 
sites are attractive {u > 0) and one expects a cluster- 
ing of sites. This will lead to a larger thermal roughness 
compared to the non-interacting case. To study the 
question whether these attractions can actually lead to 
condensation we turn to the simplified case of adhesion 
centers moving on a square lattice with Ng sites and lat- 
tice constant a. The area fraction is = N/Ns- We de- 
fine a symmetric N x N link matrix b with entries bij = 1 
if the two centers i and j are neighbors and b^j = other- 



wise. Then Xij 



number of hnks is TVl — (1/2) ^ij- 



-uibij with ui — ii(a/^||) < 1. The 
Expanding the 
excess free energy Eq. (14) to second order in ui leads to 
f3F^ « —vN]^ with interaction strength v = u\/2. This 
approximation implies only nearest neighbor interactions 
and holds for ^y < a. Hence, the statistics of the adhesion 
centers can be inferred from a two-dimensional lattice gas 
with constant density (or, equivalently, the Ising model 
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/3£b = 1.5 



/3£b = 3 




FIG. 3: Three snapshots of the lattice gas (adhesion sites are 
black squares) for area fraction <j) = 0.3 and different values of 
the binding energy e^,. The snapshots have been taken after 
10'' sweeps starting from random configurations. 



with constant magnetization) and Kawasaki diffusion dy- 
namics. This mapping constitutes the second main result 
of this Rapid Communication. However, in contrast to 
the standard lattice gas the interaction strength v is not 
constant but depends, as does , on the number of links 

We have performed Monte Carlo simulations of the 
lattice gas with Ng = 40^. In Fig. [2] we show the av- 
eraged thermal roughness as a function of the area 
fraction (j> for both the lattice gas and the ideal non- 
interacting case [Eq. (17)]. For the lattice gas the de- 



cay length is larger due to the clustering of adhesion 
centers, which implies that membrane regions can move 
further away from the substrate exhibiting a larger ther- 
mal roughness as compared to the non-interacting case. 
The decay length fluctuates with a narrow distribution 
that is Gaussian to a very good degree. For compari- 
son the thermal roughness i^j^ as obtained numerically in 
Ref. [10 is plotted. Also shown in Fig. [2] is the interac- 
tion strength vs. the area fraction. Since for a larger 
number of adhesion centers the membrane is pulled closer 
to the substrate the thermal roughness, and therefore the 
interaction strength, decrease. The effective interaction 
strength is purely entropic as one would of course ex- 
pect. Although the interaction strength is not constant 
its distribution is rather narrow and hence we compare 



its mean value (i/) to the critical value Vc ~ 1/2.27 « 0.44 
of the standard Ising model. We see that (i^) < Vc ex- 
cept for very small area fractions for which we expect the 
crossover to unconfined fluctuations. The fact that the 
simulations show no condensation (see Fig. |3] left panel) 
is therefore due to the fact that the interaction strength 
is too small. This situation changes if we allow for short- 
ranged interactions with a binding energy eb between lig- 
ands or receptors or both. Such a binding energy might 
again be entropic in nature, e.g., due to hydrophobic mis- 
match [3T] or different types of receptors with different 
lengths [H]. For large enough ly — Ui/2 + /3eb we then 
observe coarsening and condensation of bonds into large 
clusters, see Fig. [3] 

In summary, we have developed an analytic theory for 
the statistics of mobile adhesion centers (single receptor- 
ligand bonds) in adhered membranes. To this end we 
have established an upper bound to the free energy that is 
independent of details of the interactions between mem- 
brane and substrate. For sites moving on a discrete lat- 
tice the resulting minimized free energy can be mapped 
onto a lattice gas with constant density but fluctuating 
interaction strength. The mean interaction strength is 
shown to be well below the critical value of the Ising 
model. While we have focused on membrane-substrate 
interactions conceptually the same situation arises if the 
membrane is moving close to a soft surface such as 
another cell membrane. Then h{r) denotes the inter- 
membrane distance and n is to be replaced by the effec- 
tive bending rigidity k = kiK2/{ki + K2), where Ki are 
the bending rigidities of the two membranes. 
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